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Description 

This invention is concerned with measuring the rel- 
ative and absolute velocity of liquids, specifically of wa- 
ter or fluids mixed with water, as applied to oil-field op- 
erations. 

Depleted oil production zones may be rejuvenated 
by water flooding. In this process the casing opposite 
the formation under consideration is plugged off by 
packers above and below the production zone. Water 
injection tubing, nested within the casing of an injection 
well, allows injection water to flow down the tubing to 
the desired injection zone, exit the tubing through spe- 
cialized hardware, and flow in the annul us between the 
tubing and the casing. From this annulus the water en- 
ters the formation through perforations in the casing, 
pushing the oil in the formation ahead of the water-flood 
front. For good and sufficient reasons it is very important 
to the operator to monitor the velocity of the water flow 
in the tubing-casing annulus, and therefore the volumet- 
ric flow into a particular zone of the formation, through- 
out the vertical profile of the well. 

Various methods are known for measuring the wa- 
ter-flow velocity through a single conduit such as the in- 
jection tubing itself. For example, the water flow velocity 
within the injection tubing can be measured by placing 
a spinner-type flowmeter inside the tubing. However, 
that method can not, of course, be used to measure the 
flow velocity in the annulus between the tubing and the 
casing which is inaccessible to a mechanical flowmeter. 

Methods of measurement, which do not rely upon 
direct access to the flow stream in the tubing-casing an- 
nulus, have been described to measure the velocity of 
water flowing In the annulus. A radioactive tracer, for ex- 
ample, could be mixed with the injected fluid to monitor 
the progress of the fluid with detectors sensitive to radi- 
oactive decay. That method is unattractive, first because 
of the need for handling radioactive substances and, 
second, because the results are highly qualitative. 

Another method, referred to as oxygen activation, 
may be used to measure the velocity of water flow Inac- 
cessible to direct measurement. This method consists 
of continuously irradiating the oxygen nuclei of the flow- 
ing water with high-energy neutrons to generate there- 
from the unstable isotope N^® which has a half life of 
7.13 seconds. As the isotope decays, gamma rays are 
emitted which are counted by two radiation detectors 
downstream of the source that are spaced a distance, 
s. apart. The time behaviour of the recorded count rates 
as seen at the two detectors follow the relation 

=Kexp{-Xt1) (1) 

and 



Cg = K exp{-Xt2} (2) 

where C-, and C2 are the respective count rates at times 
^ t1 and t2, at distances d1 and d2 = (dl+s) from the 
source where d1 is the distance from the source to the 
near detector. K is the decay rate at time t1 = t2 = 0. and 
X is the N^^ decay constant. Using these relatbns and 
expressing time as the ratio of distance to vekx:ity, the 
velocity of activated water flowing between the two de- 
tectors, is given by 

V^=sX/(\n{0,)'\n(C^)) (3) 

Another method for using oxygen activation to 
measure the velocity of water, referred to as the impulse 
method, differs from the above method In that the source 
is turned on and off, for example In a 1 0-secondson and 
60-seconds-off pattern. In this method, a localized vol- 
ume of water in the region of the source Is activated 
when the source is on and its time of passage past the 
detector is noted as a peak in the count rate. The velocity 
of the flow maybe calculated from the time of passage 
and the known distance from the source to the detector. 

While all of the current methods that employ oxygen 
activation to measure water flow are capable of meas- 
uring the vetocity of water flow in the tubing-casing an- 
nulus. none can do so, for practical purposes, if there is 
an additional, physically separate, stream of water 
present flowing in the same direction (a co-dlrectlonal 
flow stream) at a different velocity. In that case the neu- 
tron source bombards the O"*^ nuclei of both volumes of 
water at the same time and the detectors cannot distin- 
guish between the counts from the two flow volumes be- 
cause of superposition of the decay activity. Hence, the 
measurement of the velocity of water flow in the pres- 
ence of multiple co-directional flow streams is not pos- 
sible using the current methods of oxygen activation. 

Some of the applicable patents Include US Patent 
No. 3,603.795. This patent teaches a method substan- 
tially the same as the method described by equations 
(1), (2) and (3). That patent also directs its application 
to the detection of water flow outside the casing as well 
as inside. 

A trilogy of US Patents, No. 4,032,780, No. 
4,032,778 and No. 4,035,640, teaches various aspects 
of measuring water flow in the region outside the well 
casing where water might flow between stratigraphic 
levels through channels in the casing-cement-formation 
annuli. Such channels are the result of an incomplete 
cement seal between the casing exterior and the bore- 
hole wall in the formation. Water flowing in these chan- 
nels is referred to as "behind casing water flow". 

The 780 patent teaches a method to measure the 
volume flow rate of behind casing water flow by using a 
measurement of the flow velocity and an estimate of the 
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distance R to the flow regran. The velocity of the unde- 
sired water flow is calculated from the count motes in 
substantially the same manner as described using 
equations (1), (2) and (3). 

The 778 patent teaches a relationship for the count 
rate ratio of two distinct energy regions of the gamma 
ray spectrum as a function of the distance from the gam- 
ma ray source. The distance to the flow channel is de- 
termined using this ratio. The calculation of the flow ve- 
locity is made in substantially the same manner as de- 
scribed using equations (1 ), (2) and (3). Using the meas- 
ured velocity and distance to the channel, the volumetric 
flow rate may be determined. 

The '640 patent teaches that background radiation 
due to prompt (n,7) radiation is largely avoided if the high 
energy neutron source is quickly pulsed and the meas- 
urements of activation count rates are made between 
the pulses. The linear water flow velocity of the unde- 
sired flow is calculated from the count rates in substan- 
tially the same manner as described using equations 
0), (2) and (3). 

In a paper entitled Applications of Oxygen Activa- 
tion for injection and Production Profiling in the Kuparuk 
River Field, published as paper 22130 in May. 1991 by 
the Society of Petroleum Engineers, H. D. Scott et al. 
teach use of a stationary togging instrument for meas- 
uring fluid velocity by oxygen activation using the Im- 
pulse method described above. Referring to the inter- 
ference from codirectional flows, the authors state that 
'...If flow does exist inside the tubing from zones below 
the packer, it may be difficult or impossible to quantita- 
tively interpret the data from the zone of interest be- 
cause of superposition of the flowing signals.. ." In the 
application under discussion In the Scott et al. paper, 
the annularflowvelocity to be measured is in the upward 
direction and the interfering flow is the tubing flow in the 
upward direction from zones beneath the region of the 
measurement. Later in the paper, the authors point out 
that if co-directional water flows having widely different 
velocities are present (such as 91 cm [3 feet] per minute 
and 30 m [100 feet] per minute), the two velocities may 
be individually measured by use of a third, long-spaced, 
detector for measuring the fast flow. 

None of the presently available art solves the need 
for a general and practical method capable of measuring 
the velocity of the fluid flow in the annulus between the 
inner and outer conduit in the presence of co-directional 
flow in the Inner conduit where two separate conduits 
are nested together. None of the references are directed 
to a general and practical method for measuring the flow 
velocity in the tubing-casing annulus where the differ- 
ence between the flow velocities in the tubing and the 
annulus is not great. 

According to one aspect of the invention there is 
provided a method for measuring the relative and abso- 
lute velocity of a volume of water, or fluids mixed with 
water, flowing through an annulus between an inner and 
an outer conduit pair nested in a well bore In the pres- 



ence of water flowing co-directional ly in the inner con- 
duit at the time of the measurement; a cable-suspended 
logging instrument is passed through the inner conduit, 
the instrument containing a source of high energy neu- 

5 trons, two gamma ray detectors spaced a distance s 
apart, associated signal processing and transmission 
electronics and a mechanical flowmeter. A mathemati- 
cal description is formulated for the recorded count rate 
due to the decay of the activated water flowing in the 

10 inner conduit as an analytical function of the relative ve- 
locity between the inner conduit flow and the instrument. 
The mathematical description Is formulated for each de- 
tector by the principle of velocity gauging. This analytical 
description of count rate vs. relative inner conduit veloc- 

15 ity is defined as the tubing count rate profile. 

Having determined the tubing count rate profile, the 
logging instrument is introduced through the tubing to a 
zone of interest in the well where the flow velocity of the 
water volume In the annulus between the nested Inner 

20 and outer conduits is to be measured. The logging in- 
strument is moved through the tubing in the well bore at 
a velocity v^; where O ^ v^ < v^ and where v-, is the an- 
nular flow velocity. The 0^® nuclei of the water molecules 
are activated by high energy neutrons from the source 

25 to produce the unstable isotope N""^ which decays to 
produce measurable gamma rays. The detectors meas- 
ure the total gamma ray count rate. The instantaneous 
velocity V2 of the water flow In the tubing relative to the 
instrument is measured concurrently by the mechanical 

30 flowmeter. From the tubing count rate profiles for the de- 
tectors, the tubing count rates corresponding to the 
measured instantaneous relative tubing flow velocity is 
determined and subtracted from the total gamma ray 
count rates to yield corrected count rates. The relative 

35 fluid velocity V in the annulus is determined from the ra- 
tio of the corrected count rates in the two detectors. The 
absolute velocity V^bs the annular flow is the sum of 
the cable velocity v^. and the measured relative velocity 
V. 

40 For a better understanding of the invention and to 
show how the same may be carried into effect, reference 
will now be made, by way of example, to the accompa- 
nying drawings, wherein: 

45 FIGURE 1 is a schematic diagram of an inner water- 
injection conduit, or tubing, nested inside an outer 
conduit, or casing, in a well bore and showing a ca- 
ble-suspended logging Instrument arranged for ver- 
tical movement within the inner conduit; 
so FIGURE 2 is an enlargement of that portion of FIG- 
URE 1 involving the logging instrument; 
FIGURE 3 is a plot of the measured gamma ray 
count rate as a function of relative water flow veloc- 
ity for a 5.08 cm (2 Inch) diameter channel of water 
55 adjacent to the instrument; 

FIGURE 4 is a plot of the measured gamma ray 
count rate as a function of relative water flow veloc- 
ity for an instrument centred in a 7.30 cm (2.875") 
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tubing nested in a 17.78 cm (7.0") casing; 
FIGURE 5 is a plot of the modeled count late be- 
haviour as a function of instrument velocity showing 
the outer conduit (slower velocity flow) and inner 
conduit (faster velocity flow) contributions and the 
resulting total count rate resulting from the super- 
position of the two components; and 
FIGURE 6 is a graph of the outer conduit relative 
velocity as measured using the methods of this dis- 
closure vs. the true outer conduit relative velocity 
as independently measured by a calibrated flowm- 
eter. 

Mechanical Details 

Figures 1 and 2 represent a cross section of a por- 
tion of an injection well penetrating a subsurface forma- 
tion 10 to a region associated with an injection zone 12. 
Casing 1 4 includes multiple perforations 1 6 opposite the 
porous injection zone 12. Injection tubing 18, nested in- 
side casing 14 is provided with openings 19 so that the 
injection fluid flows under pressure into the annutus 20 
between the inner conduit or tubing 1 8 and outer conduit 
or casing 14, thence into the formation via the perfora- 
tions 16 to sweep oil towards a production well (not 
shown). Packers 1 3 and 1 5 confine the injection water 
in casing 14 to a desired production zone 12. For pur- 
poses of this disclosure, the water flow volumes In the 
annulus and in the tubing are co-directional as shown 
by the arrows V1 and V2 respectively. The velocity of 
the togging instrument 22 is symbolized by arrow (ve- 
locity of cable). 

The logging instrument 22 Is a modification of a con- 
ventional neutron logging instrument such as the PDK- 
100 (Atlas Wireline Services registered service mark) 
provided by the assignee of this invention, it consists of 
an elongated mandrel 24 of suitable material, 4.23 cm 
(1.6875") in diameter, supported by a cable 26 that is 
coupled to the draw works 27 at the surface for deploy- 
ment through the inner conduit or injection tubing 18. 
The velocity of the instrument 22 as it is drawn through 
the conduit 1 8 is measured by an odometer/velocimeter 
of any well known type 29 that may be associated with 
a sheave over which the supporting cable 26 passes. 

A pulsed neutron source 26 is mounted inside one 
end of the instrument 22 and separated from the interi- 
orly-mounted near and far gamma ray detectors 30 and 
32 by a shield 34. A third gamma ray detector 36 may 
be provided. It is to be understood that the detectors 
may be mounted beneath the source as shown in Figure 
1 or above the source. The selection of the configuration 
depends upon the direction of the water flow to be meas- 
ured. A mechanical flowmeter 40 is secured to the bot- 
tom of the Instrument 22 for measuring the velocity of 
the fluid in the inner conduit relative to the instrument. 
Signal processing electronic circuitry (not shown) is in- 
stalled in compartments of the instrument to discrimi- 
nate against low level gamma ray activity in favor of the 



higher energy deriving from the activated oxygen. The 
detector count rates are digitized downhole and are tel- 
emetricaliy transmitted to the surface through suitable 
conductors In supporting cable 26 to processing and ar- 
s chival storage unit 31 at the surface. 

In operation, the neutron source is pulsed at 1 kHz 
for 28 milliseconds (ms) and is then shut off for 8 ms 
during which time the count rate measurement is made. 
The third gamma ray detector 36 may be used in con- 
junction with either of the detectors 30 or 32 to increase 
the detector spacing for the measurement. However, it 
is most often used to identify specific rock boundaries 
in the well as determined by correlation with a previous- 
ly-derived well-logging graph of natural gamma radia- 
tion. 

During a fluid velocity measurement logging run, 
four quantities are measured. The Instrument teleme- 
ters the total gamma ray count rates as measured by 
each of the two detectors 30 and 32 and the relative ve- 
locity of the fluid in the inner conduit as measured by the 
flowmeter 40, to the surface processing equipment 31. 
The logging instrument cable velocity as measured by 
the cable velocity odometer 29 is sent to the processing 
device 31 over a separate channel 33. As will be ex- 
plained now, from those four inputs, the fluid velocity in 
the outer conduit, such as the annular space 20, may 
be calculated. 

Mathematical Model 

Given the instrument located in the inner conduit 
and a volume of water or fluid mixed with water flowing 
in the inner or outer conduits, the ability to mathemati- 
cally describe the recorded count rates generated by ei- 
ther flow with an analytic function of the relative flow ve- 
locity is fundamental to the method described here. To 
examine such a description, consider the general case 
of a water stream flowing adjacent to the instrument. 
The flow stream is activated by the high energy neutrons 
emanating from the source 18. The resulting gamma 
rays from the radioactive decay of the N^® isotope are 
detected by the near and far detectors 30 and 32. The 
decay activity as measured in counts per second is de- 
scribed by 

CR = (AN)exp{-XB/v}, (4) 

where 

V = flow velocity relative to the instrument, 
CR = measured count rate corrected for back- 
ground when necessary and normalized to source 
output, 

X. = N^s decay constant. 
A, B are adjustable parameters. 

Laboratory data for near and far detectors is shown 
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in Figure 3 where the data were taken with a 5.08 cm (2 
inch) diameter flow stream adjacent to the Instrument. 
The curves 42 and 44, of the form of equation (4). have 
been constructed by optimizing parameters A and B us- 
ing weighted least squares variance minimization. The 
resulting curves provide a satisfactory description of the 
measured data. 

Laboratory data were also taken in the configuration 
where the instrument is contained in the inner conduit 
of an inner and outer nested pair of concentric conduits. 
The inner conduit was constructed of 7.3 cm (2.875") 
diameter. 2B.9N (6.5 pound) production tubing and the 
outer conduit was 17.78 cm (7') diameter production 
casing. Water could be configured to flow in the inner 
and outer conduits individually and separately or in both 
simultaneously. Figure 4 shows the measured count 
rates for the near and far detectors as a function of the 
relative velocity between the instrument and the flow. 
The data shown are for individual, separate flows in the 
inner and outer conduits. The curves show the opti- 
mized fit of equation (4) to the data. Curves 41 and 43 
apply to the near detector and curves 45 and 47 apply 
to the far detector. 

Having demonstrated that the behaviour of the 
count rate as a function of relative flow velocity can be 
described by equation (4) for flow in both inner and outer 
conduit flows, consider now the instrument deployed in 
a conduit where it is in motion co -direct ion ally with the 
water flow contained in that conduit. If the velocity of the 
instrument is the cable velocity, v^, the velocity of the 
flow in the conduit is V2^ and V2 > v^, then the velocity of 
the flow relative to the instrument is 

v = V2-v^ 

and equation (4) can be expressed as a function of v^. 
The shape of the resulting function for V2 = 6.1 m (20 
feet) per minute is shown by the curve 48 in Figure 5 
where A = 20 and B = 1 for illustration. 

In the operational configuration, an inner conduit 
containing the instrument is nested within an outer con- 
duit and there are co-directional water flows in both con- 
duits. Since the behaviour of the count rate from the out- 
er flow as a function of relative velocity can be described 
by a function of the form of equation (4), a model of the 
count rates from both inner and outer flows as a function 
of v^ is shown in Figure 5 where the individual contribu- 
tions originating from both the inner and outer flows are 
shown separately and superimposed as they would be 
recorded by the detectors. For the purpose of illustra- 
tion, the velocity of the outer-conduit flow v-, is arbitrarily 
set at 6.1 m (20 feet) per minute, that of the inner flow 
V2 is set at 12.2 m (40 feet) per minute, B = 1 for both 
flow descriptions, and the amplitude of the slower flow 
is 20% of the faster flow. The dashed curve 46 repre- 
sents the count rate contribution originating from the in- 
ner condu it. the dotted curve 45 represents that from the 



outer conduit flow, and the solid curve 50 represents the 
sum of the inner and outer contributions. Over the region 
from v^ to V2 the curves 46 and 50 are co-linear since 
the total recorded count rate is due only to the inner con- 

s duit flow. 

For v^ above 0. 203 m/s (40fpm). the instrument out- 
runs the motion of the activated fluid volumes in both 
the inner and outer conduits so no counts are recorded 
from the flows. In the range where v-| < v^ < V2, only 

10 counts from the flow in the inner conduit contribute to 
the total count rate as the instrument is outrunning the 
slower flow in the outer conduit. In the range where v^, 
< v-, < V2^ the total count rate receives contributions from 
the fluids in both conduits. 

IS The essence of this embodiment is a method for 
correcting the total count rate when operating in the 
range where v^ < v^ < V2 in order to isolate the count 
rate contribution from the outer conduit flow. This meth- 
od involves characterizing the inner conduit flow count 

20 rates by equation (4) prior to logging and subtracting the 
result from the total count rate measured during the log- 
ging run. The method of velocity gauging is used to ac- 
complish this characterization. 

25 The Tubing Count Rate Profile and the Method of 
Velocity Gauging 

The method of velocity gauging as used to analyti- 
cally characterize the gamma ray count rate from the 

30 irradiated fluid flow in the inner conduit consists of, first, 
a series of constant-cable-velocity measurements 
gauged to be contained in the cable velocity range v^ < 
Vq < Vg where the faster inner conduit flow may be meas- 
ured without interference. Then equation (4) is fitted to 

35 the resulting count rate vs. relative velocity data by 
weighted least squares variance minimization to opti- 
mize the parameters A and B thereby to formulate a 
mathematical description of the count rates as a function 
of the relative velocity between the inner conduit flow 

40 and the instrument: 

CI j = (A1 j/v)exp{-XBl /v} (5) 

45 

C2j = (A2i/v)exp{-XB2jA^} (6) 

where 

C1 j and C2j refer to the count rates due to the inner 
^ flow in the near and far detectors respectively and v is 
the relative velocity between the instrument and the in- 
ner flow as measured by the mechanical flowmeter. 
Equations (5) and (6) are defined as the tubing count 
rate profiles. Having characterized the inner conduit flow 
^ in the form of tubing count rate profiles for both the near 
and far detectors, the profiles are used to predict the 
count rate originating from the inner conduit at any ar- 
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bitrary relative velocity as measured by the flowmeter 
40 and allow the subtraction of the inner conduit contri- 
bution from the total count rate at any arbitrary relative 
velocity of Inner conduit flow and, therefore, the isolation 
of the count rates originating from the outer conduit from 
which the velocity of the annular, or outer conduit flow, 
is determined. 

On first examination it may appear that the above 
argument is circular or contradictory in that the outer 
conduit velocity, v^, is employed in the specification of 
the velocity gauging procedure, whereas the purpose of 
the measurement is the determination of v.,. The reso- 
lution is that there are some regions in each well where 
the outer conduit velocity is known to be approximately 
zero and therefore the gauging of the cable velocities 
for the purpose of characterizing the inner conduit flow 
in such regions is trivial. Such a region would be above 
the first fluid exit points in the inner tubing 18 and above 
any packers, such as 13, Figure 2, which are used to 
isolate an injection zone. It is also possible, however, to 
characterize the inner flow in a region where the outer 
flow is known to be non-zero by determining the greatest 
upper bound on the outer conduit flow velocity, v^ , based 
on calculations using the injection rates measured at the 
surface. 



characterized in the same functional form as equation 

(1): 



CI. =(D1/V)exp{-A.E1A/} 



C2^ = (D2A0exp{-XE2A/} 



(9) 



(10) 



where 



75 



V is the velocity of the outer conduit flow relative to 
the instrument, and 

D1, D2, El, E2 are laboratory-derived parameters 
characteristic of the borehole and conduit geome- 
try. 



The relative velocity V of the volume of fluid in the 
outer conduit is determined from the ratio of the correct- 
ed count rates: 



25 



V UE2-E1) 

ln(^)-ln(g|) 



(11) 



Calculation of the Outer Flow Velocity 

Following the characterization of the inner flow and 
the construction of a tubing count rate profile, stationary 
and/or continuous measurements may now be made 
over the region of interest. During these measurements, 
the flowmeter provides an Instantaneous measurement 
of the velocity of the fast flow in the inner tubing string 
18 relative to the instrument. From these measure- 
ments, a correction to the total measured count rate in 
both the near and far detectors 30 and 32 can be made: 

C1^ = Cli -Clj (7) 



C2^ = C2, - C2i (8) 

where 

the suffixes 1 and 2 denote the near and far detec- 
tors, respectively, 

CIq and C2q are the corrected outer conduit flow 
count rates, 

C1 1 and C2^ are the total measured count rates, 
and 

Clj and C2j are the Inner conduit flow count rates 
calculated from equations (2) and (3) and the flow- 
meter measurement of relative velocity. 

Since the annular flow in the outer conduit can be 



The absolute velocity V^bs of the fluid flow in the 
outer conduit is given by the sum of the measured rela- 
tive velocity and the cable velocity 

30 

Vabs-V + V, (12) 

To illustrate the validity of this method, a series of 

^ laboratory measurements of the outer conduit velocity 
were made using nested conduits consisting of 7.3 cm 
(2.875"), 28. 9N (6.5 pound) production tubing central- 
ized in 17.78 cm (7.0") casing with simultaneous co-di- 
rectional flows in both conduits. Prior to the velocity 

^ measurements, a tubing count rate profile was con- 
structed in the manner described above resulting in the 
parameters A1 , B1 , A2. and B2 describing the inner con- 
duit flow. Similarly, laboratory measurements estab- 
lished the parameters D1, El, D2, and E2 describing 

^ the outer conduit flow. Data were then taken with the 
inner conduit flow fixed at 33.5 m (110 feet) per minute 
while the co-directional outer conduit flow velocity was 
varied. The isolation of the outer conduit count rate and 
the calculation of the outer conduit flow velocity were 

^ made In the manner described above. The result is 
shown in Figure 6 where the true relative flow velocity 
in the outer conduit, as measured using a calibrated 
flowmeter on the test fixture, is plotted against the rela- 
tive flow velocity as measured using the method de- 

^ scribed herein. The vertical error bars represent the 
standard error in the measurement and are based on 
one standard deviation for equal time measurements. 



6 



11 



EP 0 645 520 B1 



12 



The results show excellent agreement between the 
measured value and the true value. 

Addrtionally. it may be observed from Figure 6 that 
the statistical precision of the measurement increases 
as the relative flow velocity decreases. The method dis- 
closed here, which employs relative measurements, 
provides the additional benefit that the precision of the 
measurement can be increased by increasing the cable 
speed to reduce the relative velocity between the cable 
and the outer flow velocity. Further, this relative meas- 
urement capability produces the significant economic 
advantage of the measurement of higher outer conduit 
flow speeds, in principle limited only by the maximum 
logging speed attainable by the mechanical logging 
equipment. 

Additionally, since the measurement of the count 
rates, cable velocity, and relative flow velocity of the in- 
ner conduit flow are instantaneous, a valid measure- 
ment of the outer conduit velocity may be made by draw- 
ing the instrument through the inner conduit at some de- 
sired velocity which may be a constant velocity, a vari- 
able velocity, or zero velocity (the instrument is station- 
ary). This capability provides the significant economic 
advantage of allowing a continuous log which allows a 
more accurate determination of the annular velocity be- 
tween perforation sets and which has not been previ- 
ously possible using traditional methods of oxygen ac- 
tivation. 

This disclosure has been written with a certain de- 
gree of specificity for purposes of illustration but not by 
way of limitation. For example, the method described 
herein is applicable also to producing welts where the 
oil is mixed with water and to any other applications 
where the measurement of the velocity of water, in the 
presence of co-directional flows, is desired. Scientists 
working in the art will conceive of variations in the meth- 
ods taught herein but which will fall within the scope of 
this invention as defined in the appended claims. 



Claims 

1 . A method for measuring the absolute velocity of a 
volume of water, or fluid mixed with water, flowing 
through the outer conduit (14) of a set of conduits 
including inner (18) and outer (14) conduits nested 
together, said velocity being measured by use of an 
instrument (22) including a source (28) of high en- 
ergy neutrons, at least two gamma ray detectors 
(30, 32) spaced a distance apart, associated signal 
processing and transmission electronics, and a me- 
chanical flowmeter (40), the method comprising the 
steps of : 

(a) formulating a tubing count rate profile for 
each said detector; 

(b) drawing said instrument through said inner 
conduit at a known velocity; 



(c) irradiating the fluid volumes flowing through 
said inner and outer conduits by high energy 
neutrons from said neutron source to produce 
the unstable N^^ isotope from the oxygen nuclei 

s contained in the water nuclei; 

(d) measuring, at each said detector, the total 
gamma ray count rates due to the decay of the 

isotope; 

(e) measuring the instantaneous velocity of the 
10 fluid flow in the inner conduit relative to the in- 
strument by said mechanical flowmeter; 

(f) defining, from the tubing count rate profile 
for each said detector, the count rates due to 
the inner flow corresponding to the measured 

IS instantaneous relative velocity between said in- 

strument and the inner conduit flow, and sub- 
tracting the so-determined count rates from the 
measured total gamma ray count rates from 
each detector to define corrected count rates; 
20 (g) determining the relative velocity of the vol- 

ume of water or fluid in the outer conduit from 
the ratio of the corrected count rates; and 
(h) adding the known velocity of the instrument 
to the measured relative velocity to obtain ab- 
25 solute velocity. 

2. A method as defined by claim 1 , wherein step (a) 
comprises the steps of: 

30 (i) passing an instrument through said inner 

conduit at each of a plurality of discreet cable 
velocities v^ that lie in a range that satisfies the 
inequality 

where v^ is the estimated greatest upper bound 
on the absolute velocity of the flow in the outer 
40 conduit and Vg is the absolute velocity of the 

water flow in the inner conduit; 

(ii) irradiating the water or fluid volumes flowing 
through said inner and outer conduits by high 
energy neutrons from said neutron source to 

45 produce the unstable N^® isotope; 

(iii) measuring the total gamnna ray decay count 
rate of the N"*® isotope at each said detector; 

(iv) measuring the velocity of the fluid flow in 
the inner conduit relative to the instrument by 

so means of the mechanical flowmeter; 

(v) formulating a mathematical descriptbn of 
the measured count rate as an analytical func- 
tion of the measured relative velocity between 
the logging instrument and the velocity of the 

ss fluid volume flowing through the inner conduit. 

3. A method as defined by claim 2. comprising: in the 
step of formulating a mathematical description, in- 
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troducing adjustable parameters in the mathemati- 
cal [(step v)] description of the measured count 
rates for each said detector; and iteratively optimiz- 
ing the adjustable parameters by weighted least 
squares variance minimization. s 

4. A method as defined by claim 1 , 2 or 3. wherein said 
tubing count rate profile predicts the count rate orig- 
inating from the inner conduit fluid flow as a function 

of the relative velocity between the logging instru- to 
ment and the inner conduit fluid flow. 

5. A method as defined by claim 1,2, 3 or 4 wherein 
said instrument is passed through said inner con- 
duit co-directionally with the fluid flow therein. 

6. A method as defined by claim 5, wherein the fluid 
flows within said nested conduits are co-directional. 

7. A method as defined by any one of the preceding 20 
claims comprising reducing the statistical uncer- 
tainty in the calculated relative outer conduit flow 
velocity by reducing the relative velocity between 

the instrument and the fluid flow in the outer conduit. 

25 

8. A method as defined by any one of the preceding 
claims wherein step (g). is defined by 

ln(^)-ln(g) 

where 

OS 

CI o and C2o = the corrected count rates, 
X = decay rate, 

D1, D2, E1. E2 are laboratory-derived con- 
stants, 

V is the fluid velocity in the outer conduit relative 40 
to the instrument. 

9. A method as defined by any one of the preceding 
claims wherein said instrument velocity is a con- 
stant velocity. 4S 

10. A method as defined by any one of claims 1 to 8, 
wherein said instrument velocity is a variable veloc- 
ity. 

so 

11. A method as defined by any one of claims 1 to 8, 
wherein said instrument velocity is zero. 

12. A method as defined by any one of the preceding 
claims when carried out with the conduits nested in ss 
a well bore and with the instrument being a logging 
instrument suspended by cable in the bore. 



Patentanspruche 

1 . Verfahren zum Messen der Absolutgeschwindigkelt 
einer Wassermenge oder eines Fluid-Wasser-Ge- 
misches, das durch die auBere Leitung (14) eines 
Rohrleitungssatzes flie3t. bei dem innere (18) und 
auQere (14) Leitungen ineinander verschachtelt 
sind, wobei die Geschwindigkeit durch den Ge- 
brauch eines Instruments (22) gemessen wird, das 
eine Quelle (28) hochenergetischer Neutronen ent- 
halt sowie mindestens zwei Gammastrahldetekto- 
ren (30, 32). die einen Abstand zueinander haben, 
zugehorige Signalverarbeitungs- und Ubertra- 
gungselektronik und einen mechanischen Durch- 
fluBmesser (40), und das Verfahren die Schritte um- 
fa3t: 

a) Formuiieren eines Steigrohr-Zahlratenprofils 
fur jeden der Detektoren; 

b) Ziehen des Instruments durch die innere Lei- 
tung mit einer bekannten Geschwindigkeit; 

c) Best rah I en der Fluidmengen, die durch die 
inneren und au3eren Leitungen flie8en, mit 
hochenergetischen Neutronen aus der Neutro- 
nenquelle, damitaus den in den Wasserkernen 
enthaltenen Sauerstoffkernen instabile N"*®- 
Isotope erzeugt werden; 

d) Messen der gesamten Gam mast rah I -Zah Ira- 
ten, die durch den Zerfall der N^^-lsotope ent- 
stehen, mit jedem der Detektoren; 

e) Messen der Momentangeschwindigkeit der 
Fluidstromung in der inneren Leitung relativ 
zum Instrument mit dem mechanischen Durch- 
flu3messer; 

f) Bestimmen - und zwar aus dem Steigrohr- 
Zahlratenprofil fur jeden der Detektoren - der 
Zahlraten aufgrund der inneren Stromung, die 
der gemessenen relativen Momentange- 
schwindigkeit zwischen dem Instrument und 
der Stromung in der inneren Leitung entspricht, 
und Subtrahieren der auf diese Weise be- 
stimmten Zahlraten von den gemessenen ge- 
samten Gammastrahl-Zahlraten aus jedem 
Detektor, um korrigierte Zahlraten zu bestim- 
men; 

g) Erfassen der Relativgeschwindigkeit der 
Wasser- oder Fluidmengen in der au3eren Lei- 
tung aus dem Verhaltnis der korrigierten Zahl- 
raten; und 

h) Addieren der bekannten Geschwindigkeit 
des Instruments zur gemessenen Relativge- 
schwindigkeit, damit man die Absolutge- 
schwindigkelt erhalt. 

2. Verfahren nach Anspruch 1 , wobei der Schritt a) die 
Schritte umfaBt: 

i) Durchfuhren eines Instruments durch die in- 
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nere Leitung jeweils mit einer Anzahl unter- 
schiedlicher Kabelgeschwindigkeiten v^, die in 
einem Bereich liegen, der die Ungleichung 

erfullt, wobei v., die geschatzte gr63te obere 
Schranke der Absolutgeschwindigkeit der Stro- 
mung in der au3eren Leitung ist, und V2 die Ab- 10 
solutgeschwindigkeit der Wasserstromung in 
der inneren Leitung; 

ii) Bestrahlen der Wasser- oder Fluidmengen. 
die durch die inneren und au3eren Leitungen 
flie3en. mit lioclienergetischen Neutronen aus is 
der Neutronenquelle, damit instabile ^-Isoto- 
pe erzeugt werden; 

ili) Messen der gesamten Gammastrahl-Zer- 
fallszahlrate der N^^-lsotope mit jedem Detek- 
tor; 20 
iv) Messen der Gescliwindigkeit der Fluidstro- 
mung in der inneren Leitung relativ zum Instru- 
ment mit Hitfe des mechanischen 
Durchflu3messers; 

V) Formulieren einer mathematischen Be- 25 
schreibung der gemessenen Zahlrate als ana- 
lytische Funktion der gemessenen Relativge- 
schwindigkeit zwischen dem Me3instrument 
und der Geschwindigkeit der Fluidmenge, die 
durch die innere Leitung flieBt. 30 

3. Verfahren nach Anspruch 2, umfassend: 

im Schritt des Formulierens einer mathemati- 
schen Beschreibung - das Einf uhren einstellba- 3S 
rer Parameter in die mathematische Beschrei- 
bung [(Schritt v)] der gemessenen Zahlraten fur 
jeden Detektor; und iteratives Optimieren der 
einstellbaren Parameter mit einer gewichteten 
Fehlerquadrat-Abweichungsminimierung. 40 

4. Verfahren nach Anspruch 1, 2 oder 3, wobel das 
Steigrohr-Zahlratenprofil die Zahlrate vorhersagt, 
die durch die Fluidstromung in der inneren Leitung 
entsteht, und zwar als Funktion der Relativge- 
schwindigkeit zwischen dem MeBlnstrument und 
der Fluidstromung in der inneren Leitung. 

5. Verfahren nach Anspruch 1,2,3 oder 4, wobei das 
Instrument in der gleichen Richtung wie die Fluid- so 
stromung in der inneren Leitung durch die innere 
Leitung gefuhrt wird. 

6. Verfahren nach Anspruch 5, wobei die Fluidstro- 
mungen innerhalb der verschachtelten Leitungen ss 
gleich gerk;htet sind. 

7. Verfahren nach irgendeinem der vorhergehenden 



AnsprOche. umfassend das Verringem der statisti- 
schen Unsicherheit in der berechneten relativen 
Stromungsgeschwindigkeit in der auBeren Leitung 
durch das Vermindern der Relativgeschwindigkeit 
zwischen dem Instrument und der Fluidstromung in 
derau3eren Leitung. 

8. Verfahren nach irgendeinem der vorhergehenden 
AnsprOche, wobei der Schritt g) bestimmt ist durch 

X{E2-E1) 

wobei gilt: 

CI0 und C2o = korrigierte Zahlraten, 
\ = Zerfa lis rate, 

D1 , D2, El und E2 sind Konstanten, die im La- 
bor bestimmt werden, 

V ist die Fluidgeschwindigkeit in der auBeren 
Leitung relativ zum Instrument. 

9. Verfahren nach irgendeinem der vorhergehenden 
AnsprOche, wobei die Instrumentengeschwindig- 
keit eine konstante Geschwindigkeit ist. 

10. Verfahren nach irgendeinem der AnsprOche 1 bis 
8, wobei die Instrumentengeschwindigkeit eine ver- 
anderliche Geschwindigkeit ist. 

11. Verfahren nach irgendeinem der AnsprOche 1 bis 
8, wobei die Instrumentengeschwindigkeit null ist. 

12. Verfahren nach irgendeinem der vorhergehenden 
AnsprOche, das ausgefOhrt wird, wobei die Leitun- 
gen in einem Bohrloch verschachtelt sind, und wo- 
bei das Instrument ein MeBinstrument ist, das an 
einem Kabel im Bohrloch hangt. 



Revendications 

1. Un proc^d^ pour mesurer la Vitesse absolue d'un 
volume d'eau ou d'un fluide m6lang6 avec de I'eau 
s'ecoulant au travers du conduit exterieur (14) d'un 
groupe de conduits comprenant des conduits int6- 
rieur (1 8) et exterieur (1 4) emboit6s conjointement, 
ladite Vitesse 6tant mesur6e par r utilisation d'un 
instrument (22) comprenant une source (28) de 
neutrons k ^nergie 61ev6e, au moins deux d^tec- 
teurs de rayons gamma (30, 32) espac6s d'une cer- 
taine distance, des dispositifs electroniques de 
transmission et de traitement de signaux associes 
et un debitmfetre m6canique (40), le proc6dd com- 
prenant les etapes consistant k : 



9 



17 



EP 0 645 520 B1 



18 



(a) formuter un profil de taux de comptage de 
tubage pour chacun desdits d^tecteurs ; 

(b) tirer ledit instrument au tfBvers dudit conduit 
interieur a une vitesse connue ; 

(c) irradier les volumes de fluide s'ecoulant au 5 
travers desdits conduits interieur et exterieur 
par des neutrons d'^nergie 6lev6e d partir de 
iadite source de neutrons pour produire I'isoto- 

pe Ni® instable d partir des noyaux d'oxyg^ne 
contenus dans les noyaux d'eau ; io 

(d) mesurer ^ chacun desdits detecteurs, les 
taux de comptage de rayons gamma totaux dus 
k la decroissance de I'isotope N^^ ; 

(e) mesurer la vitesse Instantande de T^coule- 
ment de fluide dans le conduit interieur par rap- is 
poit k I'instrument par ledit d^bitm^tre 
m^canique ; 

(f) ddfinir k partir du profil de taux de comptage 
de tubage pour chacun desdits detecteurs, les 
taux de comptage dus k I'^coulement interieur 20 
correspondant a la vitesse relative instantanee 
mesur^e entre ledit instrument et I'ecoulement 

de conduit interieur et soustraire les taux de 
comptage ainsi determines des taux de comp- 
tage de rayons gamma totaux mesur^s k partir 
de chaque detecteur pour d^finir les taux de 
comptage corriges ; 

(g) determiner la vitesse relative du volume 
d'eau ou de fluide dans le conduit exterieur k 
partir du rapport des taux de comptage so 
corriges ; et 

(h) ajouter la vitesse connue de i'instrument k 
la Vitesse relative mesurde pour obtenir la vi- 
tesse absolue. 

35 

Un precede selon la revendication 1 , dans lequel 
retape (a) comprend les Stapes consistant k : 

(1) faire passer un instrument au travers dudit 
conduit interieur k chacune d'une pluralite de 40 
vitesses de cable discretes qui s'etendent 
dans une gamme qui satisfait rinegalitS : 



ou V-, est la plus grande limite superieure 
estim6e sur la vitesse absolue de I'ecoulement 
dans le conduit exterieur et Vg est la vitesse ab- 
solue de I'ecoulement d'eau dans le conduit so 
interieur ; 

(ii) irradier les volumes d'eau ou de fluide 
s'ecoulant au travers desdits conduits interieur 
et exterieur par des neutrons k energie elevee 

k partir de Iadite source de neutrons pour pro- 55 
duire I'isotope N^® instable ; 

(iii) mesurer le taux de comptage de decrois- 
sance de rayons gamma total de I'isotope N^^ 



k chacun desdits detecteurs ; 

(iv) mesurer la vitesse de I'ecoulement de fluide 
dans le conduit interieur par rapport k I'instru- 
ment au moyen du debitmetre mecanique ; 

(v) (ormuler une description mathematique du 
taux de comptage mesure comme une fonction 
anatytique de la vitesse relative mesuree entre 
i'instrument de diagraphie et la vitesse de 
I'ecoulement de volume de fluide au travers du 
conduit interieur. 

3. Un precede selon la revendication 2. comprenant 
dans retape de formulation d'une description ma- 
thematique. I'introduction de parametres reglables 
dans la description mathematique [(etape v)] des 
taux de comptage mesures pour chacun desdits de- 
tecteurs et I'optimisation iterative des param^tres 
reglables par la minimisation de variance des moin- 
dres Carres etablie. 

4. Un precede selon I'une des revendications 1 , 2 ou 
3 dans lequel ledit profil de taux de comptage de 
tubage predit le taux de comptage provenant de 
I'ecoulement de fluide de conduit interieur en fonc- 
tion de la vitesse relative entre I'instrument de dia- 
graphie et I'ecoulement de fluide de conduit inte- 
rieur. 

5. Un precede selon I'une des revendications 1 . 2, 3 
ou 4 dans lequel ledit instrument passe au travers 
dudit conduit interieur de fa^on co-directionnelle 
avec I'ecoulement de fluide k I'interieur de ce!ui-ci. 

6. Un precede selon la revendication 5, dans lequel 
ies ecoulements de fluide dans lesdits conduits em- 
boites sent co-directionnels. 

7. Un precede selon i'une quelconque des revendica- 
tions precedentes comprenant la reduction de I'in- 
certitude statistlque dans la vitesse d'ecoulement 
du conduit exterieur relative calculee en reduisant 
la Vitesse relative entre I'instrument et I'ecoulement 
de fluide dans le conduit exterieur. 

8. Un precede selon Tune quelconque des revendica- 
tions precedentes, dans lequel retape (g) est defi- 
nie par : 

V X(E2-E1) 
ln(c^^)-ln(gl) 

OU 

CIq et C2q = ies taux de comptage corriges 

A, = taux de decroissance 

D1 , D2. El . E2 sent des constantes derivees 
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du laboratoire, 

V est la Vitesse de fluide dans le conduit ext6- 
rieur par rapport k rinstrument. 

9. Un precede selon I'une quelconque des revendica- 5 
tions precedentes. dans lequel ladite vitesse de 
rinstrument est une Vitesse constants. 

10. Un proc6d6 selon Tune quelconque des revendica- 
tions 1 k 8, dans lequel ladite vitesse de rinstrument io 
est une vitesse variable. 

11. Un proc^de selon I'une quelconque des revendica- 
tions 1 ^8, dans lequel ladite vitesse de rinstrument 

est z6ro. is 

12. Un precede selon I'une quelconque des revendica- 
tions precedentes. qui est mis en oeuvre avec les 
conduits loges dans un puits de forage at avec ('ins- 
trument qui est un instrument de diagraphie sus- 20 
pendu par un cdble dans le puits de forage. 
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